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Abstract. In the initial stage of the bottom-up picture of thermalization in 
heavy ion collisions, the gluon distribution is highly anisotropic which can 
give rise to plasma instability. This has not been taken account in the original 
paper. It is shown that in the presence of instability there are scaling solutions, 
which depend on one parameter, that match smoothly onto the late stage of 
bottom-up when thermalization takes place. 
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1. The original bottom-up picture 

In the McLerran-Venugopalan model of the color glass condensate [1], small-a; gluons 
with transverse momentum below a certain saturation scale Qs are at their maxi¬ 
mum density. When applied to a nucleus-nucleus collision at impact parameter 5, 
this scale is given by [2] 



( 1 ) 


and its value is Qs ~ 1 GeV at the relativistic heavy ion collider (RHIC). Here Ra 
is the nuclear radius, p is the nuclear number density, Nc is the number of color, a is 
the coupling and Gp is the gluon distribution of a proton. In a nuclear collision these 
gluons have a typical momentum of Qs and are freed at a time around 1/Qs after 
the initial impact. In the bottom-up picture, which is based on the observation that 
inelastic processes are no less important than elastic processes for thermalization 
[3], equilibration is driven by these hard gluons and it goes through three distinct 
stages [4]. They are a) the early times 1 < QsT < b) the intermediate times 
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a < QsT < a and c) the final stage a < QsT < a 
1.1. a) 1 < QsT < 

At the early times hard gluons dominate and because of the longitudinal expansion 
the density goes down like 

Nh - . (2) 

ayQsT) 

In the central collision region most of the gluons have small longitudinal momentum, 
P 2 <C 1, otherwise they would have wandered out of the region. But this momentum 
cannot be zero either because of broadening due to multiple scattering. Effectively 
the Pz goes through a random walk in momentum space due to the random kicks 
by other hard gluons so 

2 Ar 2 

p,, ~ A'coim^, ~ - (3) 

Pz 

where N^oi is the number of collisions a hard gluon typically has encountered at the 
time r and m\) is the screening mass square 


~ a y Sp 


,3 fhip) aNh Ql 


which effectively acts as the variance for each kick due to the much more frequent 
small angle collisions. Pz comes out to be 


Pz ~ {aNhf^^ 


(0.t)i/3 • 


Soft gluons with momentum ks are produced during these times via the Bethe- 
Heitler formula [5] to give the parametric form for Ns 


dr a(QsT)4/3 ’ 

Once produced, random scattering by other gluons energizes these soft gluons so 
that their momenta settle around ks ^ Pz- Therefore the soft gluon distribution 
becomes 

r ^ ^i_ (7) 

k^s a(Q.r)i/3- 


1.2. b) a < QsT < a 

In the intermediate times hard gluons still dominate in numbers but now fh < 1- 
This changes the scattering rate with the hard gluons so 

k‘i ~ Ncoi aQl ( 8 ) 
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is now a constant. Assuming that the screening is mainly due to the soft gluons 


2 aNs aNh 
mo -> 


Qs ’ 


one can find self-consistently that 






(9) 


( 10 ) 


1.3. c) a < QsT < a 

In the final stage most gluons are soft Ng ^ Nh. The remaining hard gluons 
will scatter with the soft gluons and lose energy via successive gluon splitting. 
Whereas in the previous stages gluon production via the Bethe-Heitler formula 
is unaffected by multiple scattering, this is no longer true as the branching gluon 
momenta now fall within the range of the Landau-Pomeranchuk-Migdal suppression 
[6]. Specifically gluon emission with momentum larger than fcnPM = m'^^/NgcattO' is 
suppressed [7]. Ngcatt is the number density of the particles that is responsible for 
most of the scatterings. In this case the formation time of the branching gluon is 
~ fcbr/fct where kt is the transverse momentum picked up by the branching gluon 
through the random kicks by the soft gluons. It can be estimated as momentum 
broadening as before but the number of collisions is now restricted by the formation 
time tf and the mean free path A, hence 

kl^m3otf/\. ( 11 ) 

The rate of branching is roughly related to the formation time via 1/tbr ~ otl'tf- 
Equating tbr with r and requiring that the soft gluon now be in a thermal bath 
Ng ~ one finds the branching momentum to be 

fcbr ~ . (12) 

Lastly equating the energy flow from the hard gluons to the soft thermal bath, the 
temperature is determined to have the linear time dependence 

Tr^a^Qlr. (13) 

We will see later on that how some of these parametric dependences are recovered 
even after instability is included into the consideration. 


2. The instability 

As mentioned previously, early on in the collision only small-a; gluons can remain 
in the central region and they have typical transverse momentum of the order of 
Qg. This describes a picture of gluons with highly anisotropic initial momentum 
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distribution. In such a situation as pointed out a long time ago [8] and more 
recently within the context of the bottom-up picture [9], it would give rise to plasma 
instability. The instability occurs because the dispersion relation for the soft gluons 
gives a negative value for the screening mass square 


m\, 


OiNh 

Qs 


(14) 


when the momentum distribution is highly anisotropic [10]. Modes with momentum 
k < rriD are unstable. For recent reviews on the topic of instability in the context 
of heavy ion collisions, one can read for example [11]. Although the growth is 
exponential in nature^ and should be very fast on the time scale of r ~ 1 /tod 
or QsT ~ 1, it is difficult for it to lead directly to equilibration because first the 
instability only produces soft particles and second Arnold and Lenaghan (the first 
paper of [10]) showed that equilibration cannot occur before QsT ~ which is 

much later than QsT ~ 1 for small 

Instability creates many soft gluons as a result. There are two possibilities for 
the system to evolve further: 

(i) When the soft particles are saturated at fs ~ 1/a further production via 
the instability will result in gluons with k ^ niD being transferred to higher 
momenta. 

(ii) Or the instability will be completely eliminated by the soft gluons at satura¬ 
tion. 


In either case, in the same spirit of the bottom-up picture, it is natural to look for 
a scaling solution which connections the end of the exponential growth due to the 
instability to final equilibration. 


3. A possible scaling solution 

The solution(s) that we propose of course still has to start with the longitudinally 
expanding initial hard gluons 


Nh 


a{QsT) 


(15) 


For gluons produced sometimes after the beginning but before r, l/Qs < tq < t 
these have [14] 


Ns{t,To) 


Ql 


ksiro) 


Qs 


a{QsT){QsToy/^ (QsTo)^^^ ’ 

(Q.ro)i/3+^/5 


afsir, To) 


(Q.t)2/3+2-5/5 


(16) 
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where Ns{t, tq) is the number density of particle produced at time tq but measured 
at r. 

For gluons produced at time r, one can write down a family of (5-parameter 
dependent scaling solutions [14] 


N, 


afs 


Ql 


Qs 


^{QsTY/ 


3-5 ’ 


(Q.t)1/3+A5 


mo 


(gsT)l/3-25/5 ’ 

Qs 

^ (Q.r)i/2-35/io 


(17) 


where 5 > 0. At 5 = 0 they coincide with the initial parametric form of the original 
bottom-up picture described in the first section. The solutions obey 


ml, 


N, 

kl 

1 


aN„ 


-{NJsY 


'-D 




with 


or 

kl 


'^col 


Nsfs- 


1 

‘T'col 




m 


D 


(18) 

(19) 

( 20 ) 
( 21 ) 


Here mn at r is determined by soft gluons produced via the Bethe-Heitler formula 
in Eq. (19). Multiple scattering ensures that these gluons gain momentum until 
they reach a value around kg given by Eq. (20). Once there they scatter once on 
the average so they are borderline as far as reaching equilibrium. 


4. The value of 6 and ml > 0? 


So far we have always given the mass mu a subscript of D which stands for the 
Debye screening mass but in all reality, we are uncertain about the sign of the 
mass square. In section 2 we pointed out that the initial momentum distribution 
was highly anisotropic, thus some soft gluon modes were unstable. Looking at 
the problem only parametrically as done in the bottom-up picture and also here 
would not help us ascertain the sign of ml- More dynamical inputs are necessary. 
One can compare the momentum distribution and from the degree of anisotropy 
deduce whether ml is negative. But the problem is more complicated than that. 
For example from Eq. (16) the contribution of the gluons produced at tq to the 
screening mass square is 


ml{T,To) 


aNs{T,To) 

fc(ro) 


g^(g«To)3V5 

QsT 


( 22 ) 


If To <C T then this contribution is clearly negative because fcs(To) is so dissimilar 
to fcs(T). However the contribution is small compared to ml, which as seen in 
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Eq. (17), has the same expression as Eq. (22) except tq is t in this case. On the 
other hand if tq ~ r then the gluon’s momentum distribution tends to be isotropic 
and it is unlikely that mjj is negative. The more difficult case is tq < t when one 
can no longer be certain when the sizes of the contributions to mjj are comparable. 
It is here that the parameter S plays a role since the ratio of the late gluon to the 
early gluon contribution goes like Larger value of 6 put more weight 

on the late-time gluons’ contribution. Better considerations and calculations are 
necessary to determine the value of <5.° 


5. Matching onto bottom-np 

The solution(s) that we proposed in Eq. (17) would not be of any value if it did not 
describe also the equilibrium phase. In fact at a time when 

Q,f ~ a-15/2(5-66) (23) 

our scaling solution becomes identical to the intermediate stage, < QsT < 

of the bottom-up picture when the basic quantities in both cases go like 

5(-l + 3<5) 3(5-3^) 

fs^a 2(5-65) ^ ~ Qs . (24) 

This is true provided 0<(5<l/3. A graphical representation of this is shown 
in Fig. 1. At this time the present solution should make a transition into the 
original bottom-up solution which remains true for the rest of the evolution as long 
as the intermediate stage of the bottom-up picture is not too affected by the initial 
presence of the instability. 

For the case when S > 1/3, we can see from Eq. (24) that fg approaching unity. 
In fact in that case at a time QsTi ~ already /s ~ 1. Much of the picture 

of the final stage of the bottom-up becomes true except that gluons produced early 
at time tq now play the part of the hard particles since Ns{t,to) > Nh and fcs(To) 
now functions as the branching momentum khv in Eq. (12) 

fc,(ro)~a4TV. (25) 


The transfer of energy is similarly via gluon branching from these gluons into the 
bath of soft gluons. Equating once again the energy flow from these gluons into the 
thermal bath with temperature T 


* ^ — 
dr dr 


Ns{t,To) 

T 


ks{To) , 


using Eq. (16) and Eq. (25) in Eq. (26) one finds 


T 


35-78(5 

Qs q;39«-io 


{Qs 


(26) 


15-36(5 

39 ( 5-10 


(27) 
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Fig. 1. Graphical representation of the matching of the scaling solution to the 
bottom-up picture at Qsf. 


At (5 = 1/3, this takes the familiar form T ~ of a linear increase of T with r 

which is characteristic of the bottom-up picture in [4]. This heating up of the bath 
of soft gluons ends when the transfer of energy to the thermal bath is complete. 
This occurs when the branching momentum fcs(To) in Eq. (25) finally reaching 
and 

~ Nhir). (28) 

At this time QsT ~ Substituting this into Eq. (27) one gets 

T ~ , (29) 

a value that is independent of 6. One sees that independent of what value <5 takes, 
as long as 5 > 1/3, the scaling solutions match up to the final stage of bottom-up 
only at the final time QsT ~ 
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Notes 


a. The speaker. 
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b. In [12] it was shown that at late times the growth changed character from an 
exponential to a linear one and in [13] for an longitudinally expanding plasma, 
the exponent was shown to be ~ ^/t as one would expect from the form of 
Eq. (4). 

c. Bddeker considered the broadening of the Pz by multiple scattering with the 

much denser unstable gluon modes instead of with the hard gluons [15]. In that 
case he found pz ~ which would suggest a value for S ~ 5/24 < 

1/3 provided that Pz takes this parametric form until the moment when the 
instability was finally eliminated. 
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